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ABSTRACT 

We propose a new dynamical picture of galactic stellar and gas spirals, based on hydrodynamic 
simulations in a 'live' stellar disk. We focus especially on spiral structures excited in a isolated galactic 
disk without a stellar bar. Using high-resolution, 3-dimensional iV-body/SPH simulations, we found 
that the spiral features of the gas in galactic disks are formed by essentially different mechanisms 
from the galactic shock in stellar density waves. The stellar spiral arms and the interstellar matter 
on average corotate in a galactic potential at any radii. Unlike the stream motions in the galactic 
shock, the interstellar matter flows into the local potential minima with irregular motions. The flows 
converge to form dense gas clouds/filaments near the bottom of the stellar spirals, whose global 
structures resemble dust-lanes seen in late-type spiral galaxies. The stellar arms are non-steady; 
they are wound and stretched by the galactic shear, and thus local densities of the arm change on 
a time scale of ~ 100 Myrs, due to bifurcating or merging with other arms. This makes the gas 
spirals associated with the stellar arms non-steady. The association of dense gas clouds are eventually 
dissolved into inter-arm regions with non-cirucular motions. Star clusters are formed from the cold, 
dense gases, whose ages are less than ~ 30 Myrs, and they are roughly associated with the background 
stellar arms without a clear spatial offset between gas spiral arms and distribution of young stars. 
Subject headings: galaxies: spiral — galaxies: ISM — ISM: kinematics and dynamics — methods: 
numerical 



L INTRODUCTION 



Spiral structures are the most prominent features 
in disk galaxies, which can be excited by tidal 
interactions with compani on galaxies (e.g . lOh et aP 
[20081: iDobbs et al.l ,2010; S truck et al.|[201lf) as w eh as 
by the central stellar ba rs (jSellwood fc Sparkd 119881 : 
iBinnev fc Tremaini 120081 (section 6.4)). Spiral struc- 
tures can also be self-excited and maintained without 
the gravitation al perturba t ions i n globally stable disks 
as proposed bv iLin fc Shul (|1964l ): They are interpreted 
as stationary density waves propagating in a stellar disk 
(see also 'Berlin & Lin 1996). In this paper, we focus 
on dynamics of the ISM in stellar spirals arisen in disk 
galaxies without external perturbations. 

In the conventional picture, the stellar density waves 
generate spiral perturbations in a galactic potential. Dy- 
namics of the interstellar medium (ISM) in the spiral 
potential have been theoretically studied since iFuiimotol 
(fr968i) . who found a standing shock, called a 'galactic 
shock', in a tightly wrapped spiral potential. Linear 
and non-linear behavior of the shock has bee n exten- 
sively studied from various perspectives (e.g.. [Roberts! 
[T96a IShu et allUlTl IWoodwardI [1971 INormanlll978ir 
The galactic spiral dust-lanes and associated star for- 
mation in the arms have since been regarded as a 
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consequence of the galactic shock. Although this is 
a widely accepted static picture of the galactic spi- 
ral, recent time-dependent, multi-dimensional simula- 
tions suggested that shocked layers in spiral potentials 
are not always dynamically stable (e.g.,lKiniet al. 200^ 
IShettv fc Ostrike3[2006[ : [Dobbs fc BoWl[2006iH Spi- 
ral shocks can be hydrodynami cally unstable ('wiggle 
instability' (Wada fc Koda 2004*)), and sometimes form 
'spurs' ((Kim fc Ostrikcr 2002). The substructures in the 
inter-arm regions are formed as a result of formation of 
clumps in the layers. In reality, various physics, such 
as self-gravity of the gas, radiative cooling, heating by 
stars, and the magnetic field could affect dynami cs and 
struc tures of the ISM in a spiral potential (see Ki m et al.l 
[20101 and references therein). Using three-dimensional 
hydrodynamic simulations, taking into account the self- 
gravity of the gas and radiative coolin g /heat ing processes 
and the supernova feedback, [Wadal (j2008[ ) found that 
the classic galactic shocks are unstable and transient, 
and that they migrate to globally quasi-steady, inhomo- 
geneous arms due to nonlinear development of gravita- 
tional, thermal, and hydrodynamical instabilities. 

One unrealistic assumption in most previous hy- 
drodynamic and magneto- hydrodynamic simulations 
of galactic spirals is, however, that the spiral po- 
tential is 'fixed'. Thus, dynamical interactions be- 

See also rece nt numerical simu lations on tidally-driven spirals 
HOh et al. 2008: Struck et al.|[20Tll 'l . 
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tween the gas and stellar components were not con- 
siderecfl. Yet, iV-body simulations have shown that 
self-induced spiral structures in stellar disks are not 
stationary (Sellwood & Carlberg 1984; Sellwood 2000; 
Sellwood & Binncv j2002t iFuchs et all 120051: iFuiii et alj 
2011; Sellwood 201^. This is also the case i n an TV- 
body disk with a central bar (jBaba et alJ 120091 ) . In any 
iV-body models, amplitudes of spiral arms change on a 
time scale of galactic rotation, or even shorter. Station- 
ary spiral waves with small pitch angles h ave ne ver been 
numerically reproduced (^Sellwood 2010^, 1201 1[ ). Each 
spiral arm in the i V-body disk is tr ansient, and re- 
currently reformed. iFuiii et all ()2011[ ) revealed a self- 
regulating mechanism that maintains multi-arm spiral 
features for at least 10 Gyrs in a pure stellar disk us- 
ing three-dimensional iV-body simulations. They showed 
that the dominant spiral modes are time-dependent, re- 
flecting a highly non-linear evolution of spiral density 
enhancements . 

What we should explore next is dynamical interactions 
between the ISM and the live (i.e. time-dependent) stel- 
lar disk. Interactions between the ISM and the stellar 
component in disk galaxies have been extensively stud- 
ied since the 1970s, and it was suggested that the ISM is 
essential fo r dynami cal structures in disk galaxies (e.g. , 
Roberts fc Shu 197a iKalnaisFT97llJog fc Solomon|[l98l 
Lin fc Bertiniil985il . However, fully non-linear evolution 
of the realistic multi-phase ISM in a live-stellar disk has 
not been investigated so far. 

In this paper, we study dynamics of the multi-phase, 
self-gravitating gas in a self-induced stellar spiral arms 
in a galactic disk using Smoothed Particle Hydrodynam- 
ics (SPH) and iV-body simulations. For comparison, a 
model using a fixed spiral potential is also presented. 
When radiative cooling and heating processes as well as 
the self-gravity of the gas are taken into account, it is es- 
sential to solve the whole three-dimensional disk; other- 
wise global non-axisymmetric, inhomogeneous structures 
and their evolution are not properly followed. We also 
implement 'star' formation from cold/dense gas and en- 
ergy feedback from supernovae, by which we can study 
what triggers formation of young stars in the ISM in spi- 
ral arms. 

One should note that the present simulations are not 
fully cosmological simulations, which are necessary to un- 
derstand a full history of galaxy formation. Therefore, 
the numerical experiments here provide us with basic 
physical insight into dynamics of spiral arms in isolated 
late-type spiral galaxies at the present universe. Dynam- 
ics of spirals developed under an effect of the stellar bar 
will appear elsewhere (Baba et al. in prep.). 

2. NUMERICAL METHODS AND MODEL SETUP 
2.1. Methods 

We used our ori ginal A^-body/hydrodyn amic simula- 
tion code ASURA ( Saitoh et al.H2008l \20m to solve the 
Newtonian equation of motions and the equations of hy- 
drodynamics using the standard SPH methods (jLucvl 
[T977t IGingold fc Monaeha^ [19771 : ISomlgel IMol) . The 

^ IDobbs fc Bonnell II2006I ) used a time-dependent grav- 
itational potential taken from an A^-body simulations of 
ISellwood fc Carlberel 119841) . and ran 3-D, two-phase (i.e., cold and 
warm gas only) Smoothed Particle Hydrodynamics simulations. 



numerical method s used here are the same as those in 
iBaba et al.] (|2009l I2010D . Here we briefiy summarize 
them. 

The self-gravity of stars and SPH particles is cal- 
culated by the Tree method with GRAPE (GRavity 
PipE), a special purpose ha rdware for A^-body simula- 
tions ()Sugimoto et al.lll990f ). Here we used a software 
emulator of GRAPE, Phantom-GRAPE (Nitadori et al. 
in preparation). The gravitational softening length is 
10 pc for both SPH and A^-body particles. 

We solve the energy equation implicitl y with the cool- 
ing fu nction (jSpaans fc Normanlll997l : IWada fc NormanI 
[200l assuming soloar metallicity, in which various cool- 
ing processed are taken into account. Heating due to 
photoelectric heating of grains and polycyclic aromatic 
hydrocarbons (PAH) is considered, assuming a uniform 
far-ultraviolet radiatio n (FUV) with th e Go = 1.0 in 
the Habing unit (Gerritsen fc Ickelll997D . We also im- 
plemented energy feedback from supernovae. 

The treatment of star formation and the heating due 
to the SN fee dback was the same as those in lSaitoh et al.l 
(|2008l |2009[) . We adopted the simple stellar popu- 
lation (SSP) a pproxima.tion, with the Salpeter initial 
mass function ([Salpeteii |1955[) and the mass range of 
0.1 — IOOM0 is assumed. The local star formation rate 
(SFR), dp^^/dt, is assumed to be proportional to the local 
gas density, Pgas, and inversely proportional to the local 
dynamical time, idyn 1/ \/Gp^: 

^ = (1) 

where C* is the dimensionless star formation effi- 
ciency parameter. The value of this parameter is usu- 
ally calibrated by the global star f ormation proper- 
ties, the Schmidt-Kennicu tt relation ([Kennicuttl 119981 : 
iMartin fc KennicuttI [20011 ). An SPH particle was re- 
placed with a star particle following the Schmidt law 
(Schmidt 1959) with a local star formati on efficiency 
C^, 0.0 3 3 in a probabilistic inanne r ([Saitoh et al.l 
[2OO8I [2OO9I: [Robertson fc Kravtsovl [20081 ) ■ if criteria (1) 
riH > 100 cm-3, (2) T < 100 K, and (3) V ■ v < 0, 
are satisffied. Note that global (galactic) star formation 
rate is not directly proportional to the local star forma- 
tion effcicncy, C* , bu t is rather controll ed by the global 
evolution of the ISM ([Saitoh et al.ll2008D . 

We implemented type-II SN feedback, where the en- 
ergy from SNe is injected to the gas around the star 
particles in the thermal energy. Each SN releases 10^^ 
ergs of thermal energy for the surrounding SPH par- 
ticles (typically 32 ± 2) particles. One gas particle 
(m ~ 4OOOM0) produces about 30 SNe during 20 Myrs. 
This energy injection produces expanding cavities of hot 
gas (T ~ 10^"'^ K) in the ISM (see also section 2 in 
[Saitoh fc Makinol[2010l ). Although the dynamics of the 
ISM on a galactic scale is mainly determined by the grav- 
itational force, local structures of the ISM on a sub-kpc 
scale are affected by the feedback from SNe. 

* (1) recombination of H, He, C, O, N, Si, and Fe; (2) colli- 
sional excitation of Hj, Ci_iv, and Oi_iv; (3) hydrogen and he- 
lium bremsstrahlung; (4) vibrational and rotational excitation of 
and H2; (5) atomic and molecular cooling due to fine-structure 
emission of C, 0"*", and O and rotational line emission of CO. 
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2.2. Live Spiral Model 

We first evolve a pure A^-body stellar disk with an 
exponential profile in a static dark matter (DM) halo 
potential whose den sity profile follows the Navarro- 
Frenk- White profile (jNavarro et al.lll997[ ). The initial 
stellar disk is ge nerated using the Hernquist's method 
(jHernq uistl 1 1 9931 ) . and then it is evolved for a few Gyrs, 
randomizing azimuthal positions of stellar particles every 
several time-step s to prevent the growth of nonaxisym- 
metric features (jMcMillan fc DehnenI 120071: iFuiii et al.l 
[Ml . Fi gure 1 shows the initial circular velocity, ve- 
locity dispersions, and Toomre's Q value as a function 
of the galactocentric distance. Once the stellar disk be- 
comes almost featureless, we added a gaseous component, 
whose total mass is 10% of the stellar disk. This is re- 
garded as an 'initial' condition in the result presented in 
§3. The model parameters of the DM halo, stellar disk, 
and gas disk are summarized in Table [TJ 

The total numbers of A^-body and SPH particles are 
initially 3 x 10^ and 10^, respectively. They are finally 
3.2 X 10^ and 0.8 x 10^, respectively. The correspond- 
ing particle masses are 11000 M© and 3200 M©, r espec - 
tively, owing to the star formation. iFuiii et ahl ()2011f ) 
showed using pure A^-body simulations that one million 
A^-body particles are necessary to reproduce long-lived 
(> 10 Gyrs) spirals, and the results do not chan g e in a 
model with A^ = 3 x 10^ particles. iSaitoh et all ()2008f ) 
explored the effect of the numerical resolutions in ASURA 
in terms of representing inhomogeneous structures of the 
multi-phase gas on a galactic scale, and concluded that 
complicated structures of the multi-phase gas can be re- 
solved with a particle mass of 3500 M©, or smaller. It 
was also confirmed that the results are converged be- 
tween models with 3500 Mq and 350 Mq. 

2.3. Rigid Spiral Model 

In order to study the effect of the live stellar disk, we 
run a test model, in which the stellar disk is replaced with 
a time-independent spiral potential. The gravitational 
potential of the stellar disk is 



Table 1 

Model parameters for each mass component (dark matter halo, 
stellar disk, and gas disk), cnfw is a concentration parameter in 
the Navarro-Frenk- White profile. 



m \ (f) A- cot is. 



(2) 



^ + zl 
1 ^ 

Kq 

on a rotating frame of the spirals. Here <I>disk.Oi "t-j ^sp, 
isp, and zq are the axisymmetric potential produced by 
the stellar disk, the number of stellar spiral arms, the 
strength of the spiral perturbation, the pitch angle, and 
the scale-height, respectively. We set igp = 30°, m = 4, 
and Esp = 0.02. The phase is i?o = 1 kpc and the pattern 
speed, fisp, is ~ 12.2 km s~-^kpc~^ {RcR = 14 kpc). 
These parameters are chosen from typical structures and 
dynamics of the live spiral model. 

3. RESULTS 

3.1. Non-steady Spiral Arms and Rs Kinematics 

Figure [2] shows four snapshots of F-band surface 
brightness reproduced by the stellar density during 975 
Myrs. Obviously spiral arms are not stationary on a time 
scale of a few 100 Myrs (see also the movie in the sup- 
plementary data). If we track one of the spiral arms, we 



Component 


Parameters 


Value 


Dark Matter Halo 


Mass 


6.3 X 10^1 Mq 




Radius 


122 kpc 




cnpw 


5.0 


Initial Stellar Disk 


Mass 


3.2 X 10^" M0 




Scale Length 


4.3 kpc 




Scale Height 


0.3 kpc 


Initial Gas Disk 


Mass 


3.2 X 10" Mq 




Scale Length 


8.6 kpc 




Scale Height 


0.1 kpc 



observe that density of the arm is no longer constant - 
a part of the arm becomes faint or denser within one ro- 
tational period. There are no steady grand-design spiral 
arms, in fact spiral arms are 'wound up'. Any arms are 
stretched by the galactic shear and attenuated. Splitting, 
bifurcating, and merging of arms are often observed. At 
the same time, new spiral arms are generated somewhere 
in the disk and as a result there are always several spi- 
ral arms in the disk. These phenomena suggest that the 
spiral arms do not propagate in the stellar disk as sta- 
tionary density waves. In the present model, the number 
of spiral arms is 4-6, but the do minant modes ar e time- 
and radially-dependent (see also lFuiii et al.l[201lL for dy- 
namical evolution of pure stellar spiral arms) . 

In the left panel of Figure [3j the rotational frequency 
of the most dominant mode of the spirals is shown as 
a function of radius at t = 1.625 Gyrs. It clearly shows 
that spiral arms mostly follow the galactic rotation fl{R) 
at any radii. They do not rotate with a constant pattern 
speed. This is always the case during the evolution, as 
shown in the right three panels, in which time evolu- 
tions of stellar densities at three different radii, i.e. R — 
4.3, 8.6, and 12.9 kpc, are plotted. The inclination of 
black/gray stripes represents rotational speed of the spi- 
ral arms, showing that they basically follow the galac- 
tic rotation (red dashed line). Each stripe changes its 
shading in several 100 Myrs, indicating local density of 
the arm does change. The radially dependent rotational 
speed of spiral arms seems to be consistent with recent 
analysis of stellar kinematics in M51 and NGC 1068, in 
which a radially declining 'pattern' speed of spiral arms 
is found using the Tremaine- Weinberg method (Merri- 
field et al. 2006; Meidt et al. 2008). The numerical 
result clearly shows that spiral arms do not propagate 
as density waves with a single pattern speed in the disk, 
but they rather behave like 'material' armf0- 

3.2. Behavior of the ISM in the Non-steady Spiral 
Potential 

In these non-stationary stellar spiral arms, the struc- 
ture and dynamics of the multi-phase ISM are quite dif- 

^ Since the arms consist of stars, and the stars are not confined 
in each arm for many rotational periods, the spiral-like density en- 
hancements also possess some similar properties in the ideal density 
waves. Detailed analysis of stellar orbits and their relationship to 
the non-steady spiral arms will be discussed elsewhere (Baba et al. 
in prep.). 
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ferent from the prevalent picture of the galactic shock. 
Figure |4] is a snapshot of cross-sections of the disk at 
three different galactic radii, revealing a typical azimuth 
distribution of stellar and densities. Relatively smooth 
curves represent azimuth distributions of the stellar den- 
sity; for example there are four arms at i? = 8.6 kpc and 
6 ~ 7 arms at i? = 12.9 kpc at t = 1.441 Gyrs. Total 
gas density, on the other hand, is more spiky. Each peak 
corresponds to a high density clouds or internal substruc- 
tures of them. It is clear that high density peaks of the 
gas tend to be associated with stellar arms. Moreover, 
this complicated structure of the gas is not steady (see 
also the movie in the supplementary data). These fea- 
tures are essentially different from the standing galactic 
shocks in a non-selfgravitating, isothermal gas in the spi- 
ral potentials. 

The concentrations of cold (T < 100 K) gas seen in Fig. 
[3] actually form spiral structures in a two-dimensional 
distribution. Figure [5] shows distribution of the cold 
gas (T < 100 K) overlaid with the stellar density. It is 
clear that gas clouds roughly trace high density regions 
of stars, and they form complicated substructures in the 
stellar spirals. The associations of dense clumps can be 
traced as a single arm connected intermittently from the 
central region to the outer disk. As clearly seen in the 
movie in the supplementary data, a grand-design spiral 
arm of the gas is temporalily formed, but it becomes 
sparse in a few '^100 Myrs. 

As shown above, if the multi-arm stellar spirals are 
developed are developed by the swing a mplification of 
density disturba nces in the galactic disk (jToomrd 119811 : 
iFuiii et all 1201 1[ ). both the gas and stellar arms follow 
the local galactic rotational velocity on average, and as a 
result the velocity of the gas relative to the stellar arms 
should be ^ 15 km or smaller. On the rotating frame 
of the spiral arms, the ISM falls into the spiral potential 
from both sides, that is there are flows of cold gas slower 
and faster than the stellar arm. The flows converge into 
condensations of cold gas near the bottom of the poten- 
tial well. The relative velocities of the gas to the stellar 
spirals are comparable to the random motion of the gas. 
In other words, the turbulent ISM 'sub-sonically' move 
near the spiral potentiaQ This is an essential difference 
in the flow in a time-independent spiral potential. 

The converging gas flows near the spiral arms are 
clearly seen in a two-dimensional map, which is distinct 
from those in a rigid spiral potential. Figure |6] compares 
the velocity field of the cold gas (T < 100 K) between 
the live and rigid spiral potential models. In the live spi- 
ral model, the velocity vectors of the gas on a galactic 
rotation frame are irregular, and they tend to converge 
into the dense gas arms. In other words, since the rela- 
tive velocity of the gas to the stellar spiral is small, the 
gas is 'trapped' by the potential troughs, forming high 
density regions. This is in contrast to what is observed in 
the rigid spiral model, where the potential is given as de- 
scribed in §2.2. The gas flows in the rigid spiral model are 

® Effective sound speed, e.g. + a'^, where Uy is veiocity 

dispersion of the gas, can be larger than tfie thermal sound speed 
Cs especially for cold gas. Therefore the clumpy gas flow can be 
regarded as sub-sonic relative to the spiral arms, in other words, 
effective pressure of the turbulent gas could dominate dynamics of 
the flow near the spiral arms, even if the local Mach number is 
larger than unity (see discussion below). 



relatively regular, and they change their directions near 
the gas arms, as typically seen in the two-dimensional 
galactic shock (jWada fc Kodall2004 ). In both the mod- 
els, supernova feedback is included, but this is not essen- 
tial to change the velocity field of the gas near the spiral 
potential, although it sometimes dissolves the associa- 
tion of gas clouds (see H'S.'d^ . One should note here that 
the gas velocities relative to the spiral potential (here 
we assume that the live stellar arms follow the galactic 
rotation) are much slower than those in the rigid spiral 
model. In the right panel of Fig. |B1 frequency distribu- 
tions of the Mach number Ai = \v\/cs in both models 
are shown, where Cg is sound velocity, Cg = {pgkTg/ ii)^/"^ 
with the gas density pg and temperature Tg of each SPH 
particle, mean molecular mass /i, and \v\ is the gas local 
velocity on a galactic rotation (live spiral model) or on a 
rotating frame of the spiral potential (rigid spiral model) . 
The sound velocity ranges from ~ 0.5 to ~ 300 km s~^ 
depending on the local gas temperature. Therefore 
represents Mach number of the ISM relative to the spiral 
potential. It is clear that the gas flows in the rigid spiral 
potential are highly supersonic (i.e. ~ 10 — 40). On 
the other hand, in the live spiral model, the median is A4 
^2. In Fig. [3 we plot the velocity field of the cold gas 
relative to the galactic rotation, i.e. mean background 
flow, in the rigid model. Although the median is A^ ~ 3 
and there are more gas with higher Mach number {A4 
> 10), the entire distributions of A4 are similar in both 
model. Therefore it is reasonable that the local density 
field of the gas in both models are similar, showing fila- 
mentary and clumpy morphology. 

If the converging gas fiows of clumpy gas collide with 
low Mach number near the spiral arms, it could cause 
weak shocks locally. Weak shocks in warm gases near the 
stellar arms can trigger forma tion of molecular clouds 
through thermal instabilities (jlnoue fc Inutsukal l200l 
iVazguez-Semadeni et al.|[201lD . 

Each stellar arm becomes faint or merges into other 
arms on a time scale of the galactic rotation, or even 
shorter. Therefore one important aspect of the gas in 
the non-steady stellar spirals is that the gas concentra- 
tions, like giant molecular associations (GMA), formed 
in the potential troughs can obtain additional kinetic en- 
ergy due to the variation of the background stellar spiral 
potential. This can be a source of dissolving the gas 
arms and generating kpc-scale non-circular motions of 
the ISM (Figure [6|). We found that most of the cold gas 
clouds have large 'epicyclic' motions of which diameter is 
comparable to typical intervals between spirals arms, i.e. 

2 — 3 kpc. In other words, the gas circulates between 
spiral arms (see the movie in the supplementary data). 

The large peculiar motions of the gas are in fact ob- 
served as proper motions of maser sources in our Galaxy 
(jReid et al.ll2009[ ). The observed pecuhar motions of the 
star forming regions are quantitatively consistent with a 
numerical barred-spiral galaxy construc ted by the sam e 
methods as those in the present paper (B aba et al.ll2009[ ). 
In Figure [H an example of such a collision between cold 
streams is shown. The sequence of snapshots during 
80 Myrs shows that concentrations of cold gas clumps, 
each of which is initially ^ 1 kpc away from the other, 
merge into one association of dense gas near the center 
of the stellar arm. The massive concentrations of cold 
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gas, which would correspond to GMA in real galaxies, 
are temporarily formed and eventually dispersed as the 
background stellar arms become faint. 

3.3. Star Formation and its Feedback in Spiral Arms 

Here we implemented 'star formation' from cold, dense 
gas. Thus, collisions of dense gas clouds shown in Fig- 
ure [6] may trigger star formation. Once the star parti- 
cles are generated from their mother clouds based on the 
star-formation criteria described in §2, they behave as 
additional collisionless iV-body particles (we call them 
'young stars') in the live stellar disk. Figure IH] shows en- 
tire structures of old (i.e. disk stars) and young stars of 
which ages are less than 30 Myrs as well as the cold gas 
at t ~ 1.6 Gyrs. The young stars tend to form clusters 
along with the spiral arms consisting of old stars. How- 
ever, there is no clear offset between young stars and the 
stellar spirals (i.e. gravitational potential), in contrast to 
the prediction of the galactic shock hypothesis. This is 
reasonable from the distribution and kinematics of cold, 
dense gas shown in Figs |4] and [5] 

In Figures. [10] and [HI evolution of a spiral arm for 
a short period (20 Myrs) is shown by density and tem- 
perature maps, respectively. The gas clouds in the white 
circle around R = 7.5 kpc are cold {Tg < 100 K), and 
are surrounded by warm (Tg ^ 10'^ — 10^ K), less dense 
gas. Around the cold clouds, several hot spots {Tg ^ 10^ 
K) caused by supernovae appear, and the dense clouds 
are partially dissolved. Figure [T2l shows a frequency dis- 
tribution of 'life time' of high density gas, for randomly 
selected 1000 SPH particles. It shows that most gas par- 
ticles do not keep nn > 100 cm~'^ for more than 10-20 
Myrs. This reflects the fact that the high density clumps 
produce massive stars, and they are eventually destroyed 
by SNe originated in the clumps. Typically a 4000 Mq 
SPH particle produces about 30 SNe during 30 Myrs. 
As seen in Figure [HI gas filaments consisted of dense gas 
complexes in the stellar spiral arms produces clusters of 
massive stars, and energy feedback from SNe originated 
in the clusters sometimes trigger expanding motions of 
the gas clouds (see also the movies in the supplementary 
data). 

3.4. Spurs 

Many spiral galaxies, especially the ones with well 
defined primary dust- lanes, ex hi bit so called spurs 
and feathers (Elmcgree^ 119801 IScoville et all 120011: 
iLa Vigne et al.l .2006, ) On the other hand, although 
the gas spiral arms shows rich substructures, we do 
not see clear periodic spurs in the inter-arm regions 
in the present model (Fig[ni). This is in contrast 
to results in previous hydrodynamic and magneto- 
hydrodynamic simulations where isothermal gas is 
evolved in a fixe d spiral potential with a constant pat- 
tern spee d, fl„ (|K im fc Ostrikeii 120021: I Wada fc Kodal 
2004; Sh cttv fc Ost rikcr 2 0061: iDobbsfc Bonncll 200^ 
Kim fc Ostrikeril2006 '). This implies that difference of ro- 
tational velocities between spiral potentials and the ISM 
is essential for forming the downstream spurs/feathers. 
In the present model, there are no strong shears or or- 
dered motions of the gas in the live spiral model, because 
both the stellar spirals and the ISM follow the galactic 
rotation, i.e. fip ~ r2(i?). If fip < il{R), as suggested 



by recent numerica l simulations of tida lly-excited spirals 
(jPobbs et al.l 120101: [Stnick et al.ll2011[ ). spurs could be 
formed. This is consistent with that spurs tend to be 
clearl y observed in well-d efined two-arm spirals, such as 
M51 (IScoville et an[200l . 

4. CONCLUSIONS 

We propose a new picture of galactic stellar and gas 
spirals, based on hydrodynamic simulations in a 'live' 
stellar disk. Here we focus on basic properties of multi- 
arm spiral structures in isolated disk galaxies without 
a central stellar bar. Tidally excited spirals in in- 
teracting galaxies are probabl y produced by the kine- 
matic density wave (IKalnaisI [1973) be yond scope of 
this paper (see e. g. lOh et al.l [2008c iDobbs et al.l 120101: 
[Struck et al.l l201l[ ).In co ntrast to the theory of stand- 
ing galactic sh ocks (e.g. iFuiimot j [1968'; Roberts '1969t 
iShu et aT]|1972[ ) in static spiral potentials caus ed by stel- 
lar density waves (|Lin fc Shdll96l IBertin fc L in 1996), 
high-resolution, three-dimen sional 7V-body/SP H simula- 
tions using the ASURA code (Sa itoh et al.ll200"8[ ) revealed 
non-steady and non-linear evolution of gas and stellar 
spirals. Morphology of the spiral structures looks simi- 
lar to that seen in late-type spiral galaxies. Important 
features found in the numerical simulations and their im- 
plications are summarized as follows: 

1. Rotational speeds of stellar spiral arms follow the 
galactic rotation. Grand-design spirals therefore 
cannot maintain their shapes over a rotational pe- 
riod. As a result they are wound up and torn off 
by the galactic shear. Spiral arms are recurrently 
formed, and often merge into other arms, by which 
local densities of the arm are no longer constant. 
This non-steady nature of spiral arms has been 
reported in 2-D and 3-D A^-body s imulati ons of 
pure stellar disks (e.g.. iSellwood fc~Carlbera 119841 : 
ISellwoodI [20Tlt IFuiu et al.l 120111). We confirmed 
that this is also the case in a stellar disk with the 
ISM. 

2. The gas component on average also follows the 
galactic rotation. As a consequence, the Mach 
number of the gas relative to spiral potentials is 
around '~ 2, while on the other hand, it is more 
than 10 in the rigid spiral model. Flows of cold, 
dense gas converge toward the stellar arms, result- 
ing in formation of massive associations of the cold 
gas and star formation. This is an essential differ- 
ence from the classical galactic shock, where the gas 
passes through the potential and forms a oblique, 
global shock along the spiral potential. 

3. Most of the cold and dense gas clouds are associ- 
ated with the time-dependent stellar spirals. The 
morphology of the gas spirals look similar to ob- 
served dust-lanes in multi-arm spiral galaxies. We, 
therefore, suspect that dust-lanes in late-type spi- 
ral galaxies, especially those in multi-arm spirals, 
are not galactic shocks. 

4. Dense gas regions in the bottom of stellar arms are 
eventually dissolved into inter-arm regions. They 
are affected by variation of the background stel- 
lar spiral potential on a time scale of ~ 100 Myrs. 
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This can be a source of dissolving the gas arms and 
generating non-circular motion of the ISM. Energy 
feedback from supernovae also contributes to dis- 
solve the complex of dense gas clumps. The entire 
motions of the cold gas on a galactic rotating frame 
are 'epicyclic', but whose radii are comparable to 
intervals between spiral arms. 

5. The mechanism presented here may also explain 
complicated sub-structures of dust-lanes and dis- 
tribution of_jnokculBr_gasjii spiral galaxies, suc h 
as M33 (jOnodera et a l. 2010: Graticr et J]|20100, 
or our Galaxy. In fact, Baba ct al. (201^ showed 
using the same method that features in observed 
position-velocity diagrams of CO and HI in our 
Galaxy are identified as clumpy gas spiral arms. 

6. In the model, stars are formed from cold, dense 
gases. Consequently, young stars, whose ages are 
less than 30 Myrs, form clusters, and they are as- 
sociated with the background stellar arms. We do 
not see clear spatial offset between gas spiral arms 
and distribution of young stars. 

The results above imply that it is essential to solve both 
live stellar disk and the multi-phase ISM self-consistently 
in order to understand dynamical structures in spiral 
galaxies and their evolution. If we change parameters of 
model galaxies, such as the rotation curve and disk mass, 
we can reproduce spiral galaxies with different morphol- 
ogy. It is also important to realize that structures in 
spiral galaxies should be 'time-dependent'. Some differ- 
ences in morphologies of isolated spiral galaxies, such as 
the number of spiral arms and substructures in spirals 
could be explained by evolutional features. Barred spi- 
rals, in which m = 2 arms are driven by a stellar bar, will 
be separately discussed elsewhere (Baba et al., in prep). 
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R [kpc] R [kpc] R [kpc] 

Figure 1. Initial condition of the stellar disk models in the live spiral model, (left) The circular velocity and contributions of individual 
component as a function of the galactocentric distance, R. (center) The velocity dispersions in the cylindrical coordinates {R, (j>, z) as a 
function of R. (right) Toomre's Q value as a function of R. 
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Figure 2. Evolution of stellar spirals on a rotation frame at 
-R = 8.6 kpc. Gray-scale represents V-band surface luminosity cal- 
culated f rom stella r density with a population synthesis model, PE- 
GASE.2 l IFioc fc Rocca-Volinerange 1999), assuming a solar metal- 
licity and the Salpeter initial mass function. 
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Figure 3. (left) Rotational frequency of the most dominant mode (m = 4) of the spirals is shown as a function radius at a snapshot 
of t = 1.625 Gyr. The gray scale represents the relative amplitude of the mode. The solid curve represents the circular frequency 
(right three panels) Kinematics of stellar arms at R = 4.3, 8.6 and 12.9 kpc. Black/gray stripes represent how the azimuthal positions of 
high density regions of stars move. The red solid lines represent the galactic rotation at each radius. 
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Figure 4. Azimuth distributions of stellar (red) and total gas (blue) density contrasts at i? = 4.3,8.6 and 12.9 kpc at t = 1.441 Gyr. See 
also the movie in the supplementary data for time evolution, where azimuthal density distributions are plotted on a rotating frame of a 
galactic rotation at each radius. 
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Figure 5. Cold gas (T < 100 K), whose surface density is rep- 
resented by the eolor, is overlaid on a stellar density shown by a 
gray-scale. See also the movie in the supplementary data. 
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Figure 6. (loft) Kinematics and surface density of tlie cold gas (T < 100 K) in the live spiral model. Arrows represent velocities of tfie 
gas relative to the galactic rotation, therefore they show approximately relative velocities to the local stellar spiral arms (see Fig. [Sj. Color 
represents gas density, (middle) Same as the left panel, but for the rigid spiral model. The velocities are on a rotating frame of the spiral 
potential. Note that unit length of the vectors in the two panels differ, and the velocities are much faster in the right two models than 
those in the live model, (right) Frequency distributions of the Mach number of SPH particles A4 = \v\/cb (see the text for the definition) 
in the live spiral model (black) and rigid spiral model (red). 
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Figure 7. Same as the middle and right panels of Fig. |6] but the velocities are defined on the local galactic rotation. 
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Figure 9. (Left) V-band surface luminosity calculated from stellar density with a population synthesis model. (Right) Cold gas (T < 100 
K) density and star particles, whose ages are < 30 Myrs represented by light blue color, formed from cold, dense gas are overlaid on the 
left panel. 
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Figure 11. Same as Fig. 1101 but temperature distribution. The hot regions (red) are generated by supernova explosions, (see also the 
movie in the supplementary 
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Figure 12. Frequency distribution of 'life time' to keep number densities, njj > 10, 10^'^ and lO'^ cm~^. 



